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Abstract 

The master equation describing the completely positive time evolution 
of a uniformly accelerated two-level system in weak interaction with a 
scalar field in the Minkowski vacuum is derived and explicitly solved. 
It encodes the well known Unruh effect, leading to a purely thermal 
equilibrium state. This asymptotic state turns out to be entangled when 
the uniformly accelerating system is composed by two, independent two- 
level atoms. 
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1. INTRODUCTION 



When a particle detector moves with a uniform acceleration through an external field 
in its vacuum state, spontaneous excitation can occur; indeed, the detector behaves as if 
it were in a thermal bath, with temperature (the "Unruh temperature") proportional to 
its proper acceleration. This phenomenon is known as the "Unruh effect" . [1-6] Although 
originally presented for detectors in interaction with massless, free scalar fields, the effect 
has been later confirmed in the case of massive and higher spin fields, in arbitrary space- 
time dimensions (with appropriate gray-body corrections to the thermal bath spectrum) , [5] 
and even for arbitrary interacting fields, using the algebraic formulation of quantum field 
theory. [7, 8] Experimental verification of the phenomena is hard to be realized in prac- 
tice, since measurable thermal effects require very high accelerations; nevertheless, clever 
proposals involving circular trajectories have been discussed and might get realized in the 
future. [9] 

Most treatments of the Unruh effect focus on the study of the spontaneous excitation 
of an accelerated "DeWitt detector", [10] a non-relativistic n-level system, linearly coupled 
to the external relativistic field. What is then studied is the excitation rate, i.e. the 
probability per unit time of a spontaneous transition from the ground state to one of its 
excited states. This probability clearly vanish for a detector at rest, but turns out to be 
non-zero, and proportional to a Planckian factor, for a uniformly accelerated one. The 
same physical result can be obtained by examining the vacuum state of the external field 
as seen by the accelerating detector: it turns out to be thermal, and described by the same 
Planckian factor. 

This result, known as the "thermalization theorem", [5] implies that the accelerat- 
ing detector behaves like an open system, i.e. a system immersed in an external heat 
bath. [11-15] In this respect, the study of its full dynamics and not only of its asymptotic 
excitation rate is of great physical interest. As for any open system, it can be obtained 
from the complete time-evolution describing the total system (detector + external field) 
by integrating over the field degrees of freedom, that in fact are never observed. 

In the following, we shall explicitly derive and study such reduced dynamics, restricting 
our attention to the case of a two-level detector in interaction with free, massless, scalar 
fields. As in the standard analysis of the Unruh effect, this simplified situation is not 
really restrictive since it is able to capture all the main features of the thermalization 
phenomena, without the algebraic complications needed to treat the case of higher-spin or 
massive fields. 

In addition, we shall assume the interaction between the accelerating two-level system 
and the scalar fields to be weak; [12-14] while encompassing the most common physical 
situations, this hypothesis allows a rigorous, mathematically sound, derivation of the evo- 
lution equation describing the dynamics of the moving detector. [16, 17] The corresponding 
finite-time evolution takes the form of a one-parameter semigroup of completely positive 
maps, [11-14] whose action on the state of the system can be explicitly obtained in closed 
form. Explicit expression for the probability transitions among the states of the moving 
detector can then be obtained: they contain appropriate Planckian factors, so that the 
known forms of the spontaneous excitations rates are reproduced. Indeed, one finds that 
the equilibrium configuration of the detector is exactly thermal, with temperature equal 
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to the Unruh temperature. This is just another manifestation of the "thermalization the- 
orem": a single, two-level system accelerating through a vacuum field is asymptotically 
driven to a thermal equilibrium with the bath. 

A physically much more interesting asymptotic situation can be obtained when the 
uniformly accelerating system is formed by two, independent two-level systems, in weak 
interaction with the same scalar field. The corresponding master equation that describes 
the reduced dynamics of the two detectors can be explicitly written generalizing the results 
obtained in the case of a single detector. It again generates a semigroup of completely 
positive maps, whose asymptotic equilibrium state can be easily determined. In general, 
it turns out to be an entangled state, even in the case of a separable initial state. 

It is known that, in certain circumstances, heat baths belonging to a specific class 
can enhance entanglement rather than destroying it: [18-23] two, mutually non-interacting 
systems immersed in one of these baths can then become quantum correlated; this can 
happen both at finite time and in the asymptotic regime. It is remarkable that the thermal 
bath seen by uniformly accelerating systems precisely belongs to the mentioned class. In 
turn, this result may suggest a new possibility for an experimental test of the Unruh effect: 
use appropriate quantum optics devices to detect the asymptotic entanglement generated 
by the uniform acceleration. 

The next Section will be devoted to the discussion of the weak coupling limit and 
the derivation of the master equation describing the dynamics of the accelerating system 
in its comoving frame. The finite time evolution it generates will be analyzed in detail 
in Section 3 for a single two-level system. The case of a detector composed by two, in- 
dependent two-level systems will be instead treated in Section 4, and the possibility of 
generation of mutual quantum correlations analyzed; in particular, the entanglement con- 
tent of the asymptotic equilibrium state will be determined by computing its corresponding 
concurrence. Finally, Section 5 contains further physical considerations on the open sys- 
tem approach to the Unruh effect, while the Appendix illustrates some technical issue 
connected with the derivation of the master equation. 



2. MASTER EQUATION 

As explained in the introductory remarks, we shall use well-known techniques de- 
veloped in the study of open quantum systems to analyze the behaviour of an idealized 
detector immersed in a scalar vacuum field, moving along a prescribed, uniformly acceler- 
ating trajectory. 

A frame-independent description of the dynamics of such a system would require a 
covariant formulation of the evolution equation of quantum mechanics. ^ Although rather 
cumbersome in practice, such a description exists: it has been introduced in the early days 

t Combining quantum mecahnics and relativity in a fully consistent way is notoriously 
problematic; for recent reviews on different aspects of this question, see Refs.[24, 25]. 
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of quantum field theory, and it is based on the so-called Schwinger-Tomonaga equation. [26] 
It extends the standard, hamiltonian evolution to a local functional dynamics, where the 
time variable is replaced by a collection of three-dimensional spacelike hyper surfaces, upon 
which spacetime has been sliced. [15, 27] 

As an alternative to this rather complicated approach, one can fix ab initio a reference 
frame and study the dynamics of the system with respect to the natural (albeit in gen- 
eral only locally valid) time variable there defined. This choice is particularly convenient 
when a global time-variable is available, as in an inertial or in an uniformly accelerated 
frame, [3] since it allows drastic simplifications in the analysis and interpretation of the 
system dynamics. The price to pay is the lack of manifest covariance and thus of a simple 
way of translating a physical description obtained in one reference frame to that derived 
in a different one. 

Both approaches are in line of principle always available, although, for practical rea- 
sons, one usually works in the "laboratory" frame. However, in certain physical situations, 
the use of the covariant formulation may in practice be precluded. This is indeed the case 
of a system in interaction with an external environment: in fact, a consistent, reduced 
dynamics can be derived only within reference frames for which environment correlations 
decay sufficiently fast. [12-15] Notice that, in most cases, this is not really a serious draw- 
back, since these preferred reference frames are precisely the ones used in the experiments. 

For the system under study, since one is interested in the behaviour of the accelerating 
detector, the natural reference frame to be adopted is that comoving with it: its dynamics 
will then be described as an evolution in the proper time t. In this reference frame, the 
detector is always at rest, and without loss of generality, it can be positioned at the origin 
of the comoving spatial coordinates. 

In the inertial reference frame, with Minkowski coordinates (x° X X X 3 ), the detec- 
tor is instead seen following an hyperbolic trajectory. More specifically, in order to better 
conform with the standard open system paradigm adopted below, we shall assume that the 
composed system {detector + external fields} be initially prepared in a factorized state, 
with the detector at rest and the fields in their vacuum state. The detector starts moving 
at time t = x° = 0, after which it follows the path: 

x°(t) = — sinhat , 
a 

x 1 (t) = - cosh at , ( 2-1 ) 
a 

x 2 (t) = x 3 (t) = , 

with a the constant proper acceleration. ^ Our aim is to derive and study the appropriate 
reduced "master equation" that describes the quantum evolution of the detector following 
this trajectory for positive t. 

To this aim, the details of the detector internal dynamics result irrelevant: one can 
then choose to model it as a simple two-level atom, i.e. as a non-relativistic, quantum 

t This situation looks more physical than the usually adopted one for which the detector 
follows the hyperbolic motion for all times. 
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mechanical system, that can be fully described in terms of a two-dimensional Hilbert 
space. With respect to a fixed, arbitrary basis in this space, states of the system will be 
represented by a 2 x 2 density matrix p, i.e. an hermitian, pt = p, normalized, Tr[p] = 1, 
operator, with non-negative eigenvalues, det[p] > 0. 

In absence of any interaction with the external scalar fields, the atom internal dynam- 
ics will be driven by a 2 x 2 hamiltonian matrix H$, that in the chosen basis can be taken 
to assume the most general form: 

3 

#s = |5>^ = f n-*, (2.2) 

i=i 

where cx^, i = 1,2,3 are the Pauli matrices, n^, i = 1,2,3 are the components of a unit 
vector, while u> represents the gap between the two energy eigenvalues. 

As mentioned in the Introduction, the interaction of the atom with the external scalar 
fields is assumed to be weak; it can then be described by an hamiltonian H' that is linear 
in both atom and field variables: 

3 

If' = J> M , (2.3) 

with oo the 2x2 unit matrix. As explicitly indicated, the interaction is effective only along 
the trajectory (2.1), as the atom is assumed to be an idealized point, without size.t The 
operators (or better, operator valued distributions) $ M (x) represent the external fields and 
satisfy the massless Klein-Gordon equation. They can be expanded as: 

N 



M*) = E [x a A~ ] ( x ) + (x a S<P^(x)} , (2.4) 



a=l 

in terms of positive- a + \x) and negative-energy (f)^\x) field operators relative to a 
set of N independent, massless, free scalar fields, with total hamiltonian H®; the complex 
coefficients "embed" the field modes into the two-dimensional detector Hilbert space and 
play the role of (generalized) coupling constants. The explicit form of the field hamiltonian 
H$ need not be specified; it suffices to know that once transformed in the inertial reference 
frame it leads to the standard expansion of </>i + ' ) , </>i in terms of Minkowski creation, 
annihilation operators, respectively. 

In the comoving frame, the total hamiltonian H for the complete system {atom + 
external fields} can thus be written as 

H = H s + H' + H$ . (2.5) 



t For a detailed discussion on the physical motivations justifying this simplifying as- 
sumption, see Ref.[5]. 



5 



It generates the evolution in t of the corresponding total density matrix p tot , 



d -^ = -iL H [pUt)], (2-6) 

starting at t = from the initial configuration: ptot(O) = p(0) ® |0)(0|, where |0) is the 
Minkowski field vacuum state, as seen in the comoving frame; the symbol Lh in (2.6) 
represents the Liouville operator corresponding to if, 

L H [.] = [H,.]. (2.7) 

The dynamics of the atom is then obtained by summing over the field $ degrees of freedom, 
i.e. by applying to ptot(^) the trace projection operator P: 

pit) = P[ Ptot {t)] =Tr $ [p tot (t)] =Tr* [e"^ (p(0) ® |0> (0|) e ltH ] ; (2.8) 

it is generated by an equation that is the result of the action of P, and its complement 
operator Q = 1 — P, on both sides of (2.6). It can be conveniently written as an integro- 
differential equation: [12] 



dp(t) 
dt 



-iL Hs [p(t)}+ ( dsK, 8 [p{t-s)] , (2.9) 
Jo 



where the time-dependent kernel JC S is defined as a formal power series expansion in the 
interaction hamiltonian H' , and it is explicitly given by the expression: 

lC s [p} = -Tr${L H > Q e - isQL " Q Q L H ,[p ®|0)<0|]} . (2.10) 

The resulting finite evolution map p(0) h- > p(t) is in general very complicated, devel- 
oping irreversibility and memory effects; indeed, the atom state p(t) at time t as given by 
(2.9) depends not only on the initial state p(0), but also on all states p(s), with s < t. This 
is a general result, valid for any reduced dynamics: nevertheless, the form of the master 
equation (2.9) can be further simplified on the basis of additional physical considerations, 
that in the present case amounts to the requirement of a weak coupling between the moving 
atom and the external fields. [11-15] 

On general grounds, one expects the memory effects in (2.9) to be negligible when the 
ratio r/r$ between the typical variation time r of p(t) and the decay time r$ of the field 
correlations is large. This is precisely the result of the "weak coupling limit" procedure: [13] 
it first amounts to a rescaling of the interaction hamiltonian H' by a dimensionless coupling 
constant g, so that the first contribution in the r.h.s. of equation (2.9), representing the 
unperturbed motion of the atom, is of order one, while the remaining piece, taking care 
of the interaction with the external field, becomes of order g 2 . This implies that the 
evolution of the state p(t) develops on time scales of order 1/g 2 ; in order to obtain a 
consistent physical description, one then needs to appropriately rescale the time variable, 
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t -> t/g 2 . After this, g can be safely taken to be vanishingly small, and the evolution 
equation attains a well defined limit. 

The procedure that we have briefly outlined can be given a precise mathematical 
meaning: [16, 17] provided the field correlations decay sufficiently fast at large time sepa- 
rations, one can rigorously prove that in the weak coupling limit the formal expression in 
(2.9), (2.10) converges to a differential equation for the reduced density matrix p(t), that 
is local in time. It takes the so-called Kossakowski-Lindblad form (see the Appendix for 
details): [28, 29] 

^ = -i[H eS ,p(t)]+jC[p(t)) , (2.11) 

with 

1 3 

£[p] = ^ £ aij[2<TjP<Ti - <Ti<Tj p - p<Ti<Tj] . (2.12) 

i,J=l 

The corresponding finite-time evolution maps ^ t generated by this equation, p(0) h- > pit) = 
7t[p(0)], form a one-parameter semigroup of trace-preserving transformations: r )t°ls = 

lt+s, t,8>0. 

The effective hamiltonian H e Q and the coefficients of the 3x3 Kossakowski matrix 
depend on the Fourier transform of the field vacuum correlations (Wightman functions), 
evaluated along the trajectory (2.1): 



aw(A) 



/oo 
dte iXt (0\$^(x(t))^(0)\0) . (2.13) 
-oo 



Being invariant, these correlation functions can be computed in any reference frame. Re- 
calling that the field variables (j)^ in (2.4) are all independent, one finds: 



N 



m,(x)^(y)\o) = J2^{x:Y G(x-y) , (2.14) 



a=l 



where G(x — y) is the standard four-dimensional Wightman function for a single scalar 
field, that, with the proper is prescription, can be written as^ 

G(x) = [ 9(k°) d(k 2 ) e^- £k ° = ~ fQ .] 2 7^0 • (2-15) 

J (2tt) 3 4tt 2 (x° - is) 2 - (x) 2 

Its Fourier transform along the trajectory (2.1) can be easily evaluated through a contour 
integral: [5] 

/oo 1 A 

Jte^G(x(t)) = - T -^, (2.16) 



^ The existence of the weak-coupling limit is guaranteed by the convergence of the 
integral J °° dt \G(x(t))\ (1 + 1) 5 , for some 5 > (for details, see Ref.[16]); this is assured 
at infinity by the exponential fall off of G(x(t)), and at zero by the is prescription. [5] 
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where flu = l/TJy, with T v = a/27r, the so-called Unruh temperature. With the help of 
the 3x3 hermitian matrices (see Appendix), 



V4° = rii 



n 



j i 



^iP = 2^ ij ~ niUj ± ie ^ kUk ) 



(2.17) 



it proves convenient to define the transformed coupling coefficients x^p a = J^jXj^jfj 
with £ = 0,+,—, together with their corresponding complex conjugate ones xP a = 
^7 {XjY ^ji^'i ^ means of them, the Kossakowski matrix in (2.12) explicitly reads: 



N 

= E [$(°) xf )a xf a + X, (+)a ^ +)a + X, (_)a ^" )a ] • (2-18) 



Being the sum of three manifestly positive terms, the hermitian matrix ciij turns out to 
be positive. As a consequence, the one parameter family of transformations 7t generated 
by the equation (2.11), (2.12) is composed by completely positive maps. [12-14, 30] As 
well known, this property assures the positivity of the evolved density matrix p(t) in 
any physical situation, thus guaranteeing the correct interpretation of its eigenvalues as 
probabilities. (For discussions on this relevant point, see Refs.[31-34].) 

This result is non-trivial and shows the importance of adopting a physically consis- 
tent and mathematically precise procedure, the Davies weak coupling limit, in deriving 
the reduced evolution equation. [35] Indeed, note that direct use of the standard second or- 
der perturbative approximation in the original master equation (2.9) (as adopted in [36]) 
produces a finite time evolution for p(t) that in general does not preserve the positivity of 
probabilities. 

Besides producing the non-unitary evolution term (2.12), the coupling with the ex- 
ternal fields 3?^ induces also a correction to the system hamiltonian (2.2), the so-called 
Lamb shift Hu- the complete hamiltonian is now if e ff = Hs + Hu- As for the Kossakowski 
matrix, this additional shift Hl can be expressed in terms of the field correlations along 
the accelerating trajectory. Introducing together with the Fourier transforms (2.13) also 
their corresponding Hilbert transforms: [12, 14] 



P 



in 



fl^u(X) 

where P denotes principal value, one finds: 



dz ^ 



A 



(2.19) 



(2.20) 



with 



3 = 1 
3 

+ E tikl 
k,l,r,s=l 



«,o(0)- fl j(0) -fl jO (0) 



(0) 



(2.21) 



(-),/,(+) 
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This expression is however formal and requires renormalization. It involves the following 
integral transform of the scalar Wightman function (compare with (2.19)), 



P f°° 
K A = — / dz 

™ J-oo 



G{z) 
z-X 



(2.22) 



which, recalling (2.16), can be split as 



K{\) = 



p r 

2^iJ 



dz 



+ 



P 



z-X 2tt 2 



dz 



1 - e^ 2 



z + X z — X 



(2.23) 



into an inertial and an acceleration dependent piece. Although not expressible in terms of 
elementary functions, [37] the acceleration dependent second term is a finite, odd function 
of A, vanishing as j3u becomes large, i.e. for a vanishing acceleration. The first contribution 
in (2.23) is however linearly divergent. As a consequence, despite some cancellations that 
occur in (2.21) (see below), the Lamb contribution Hl turns out to be infinite, and its 
definition requires the introduction of a suitable cutoff and a renormalization procedure. 

This is a well known fact, and has nothing to do with the weak-coupling assumptions 
used in deriving the evolution equation (2.11), nor with the specific situation of an accel- 
erating atom: as the splitting in (2.23) shows, the Lamb shift would be infinite even for 
an atom at rest. Rather, the appearance of the divergences is due to the non-relativistic 
treatment of the moving two-level atom, while any sensible calculation of energy shifts 
would have required the use of quantum field theory techniques. [38, 39] 

In our quantum mechanical setting, the procedure needed to make the Lamb con- 
tribution Hl well defined is therefore clear: perform a suitable acceleration independent 
subtraction, so that the expression in (2.21) reproduces the correct quantum field theory 
result when the atom is at rest. However, since we are interested in analyzing the effects 
due to the uniformly accelerated motion of the atom, we do not need to do this explicitly. 
In the following we shall therefore ignore standard, acceleration-independent hamiltonian 
contributions in the evolution equation (2.11) and concentrate on the phenomena induced 
by the motion of the atom. 



3. SINGLE ACCELERATING ATOM AND DECOHERENCE 

We shall now explicitly discuss the dynamics of the accelerated atom as described by 
the evolution equation (2.11), (2.12). In order to simplify a bit the treatment and be able 
to write explicit expressions for the evolved state p{t), we shall assume that the coupling 
coefficients introduced in (2.4) satisfy the further condition: 

N 

: (3.1) 

a=l 
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in this way, the field correlations in (2.14) become diagonal; in the following, for sake of 
simplicity, we shall set to one the proportionality coupling coefficient. Then, the sum in 
(2.18) can be explicitly performed and the Kossakowski matrix takes the general form 



Q>ij — A Sij 



where 



B 

C 



g(u) + g(-u) 
g(oj)-g(-uj) 
2g(o)-g(u)-g(-u) 



iBeijkUk + C riiUj , 
-l + e -Pu«>- 



(3.2) 



ui 

47T 

UJ 

4n 



1 



-/3 V U> 



(3.3) 



UJ 

4tt 



1 + e 



/3uU> 



Similarly, also the Lamb shift contribution b{ in (2.21) simplifies, and become directed 
along the unit vector ft; the effective hamiltonian in (2.11) can then be written as in (2.2), 



H. 



eff 



n ■ a 



in terms of a renormalized frequency 

tt = uj + i[JC(-u) - 1C(uj)] 



(3.4) 



(3.5) 



As explained at the end of the previous section, a suitable acceleration independent subtrac- 
tion has been implicitly included in the definition of the combination JC(-uj) — K,(oj), which 
otherwise would have been logarithmically divergent (compare with the result (2.23)). 

In order to discuss the properties of the solutions of (2.11), (2.12) it is convenient to 
express the density matrix p in terms of the Pauli matrices; recalling the normalization 
condition Tr[p] = 1, one has the standard expansion: 



1 

P = g i 1 + XI Pi ^) • 

i=l 



(3.6) 



Then, the evolution equation (2.11) can be conveniently rewritten as a Schrodinger-like 
equation for the coherence (Bloch) vector \p(t)) of components {pi(t), p 2 (t) 7 P3(t)}:[14] 



§ i \p(t)) = -2H\p(t)) + \ V ). 



(3.7) 



The constant vector \rf), with components rji = —ABrii, i = 1,2, 3, comes from the imagi- 
nary part of the Kossakowski matrix (3.2), while the 3x3 matrix H includes contributions 
both from H e ff and the real part of a^ : 



n 



a b + 3 c — 0,2 
b - 3 a p + flt 
c + 2 P - fti 7 



(3.8) 
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(3.9) 



where Oj = (0/2) n^, i = 1, 2, 3, and 

a = 2A + C{n\ + nf) 
a = 2A + C(n 2 1 + ttJ) 
7 = 2A + C{n\ +n\) 

The matrix TL is non singular; indeed, its eigenvalues can be explicitly determined: 
Ai = 2A, A± = (2 A + C) ± in/2. Further, their real parts are positive, so that for large 
times \p{t)) reaches an equilibrium state |poo).[40] This asymptotic state can be easily 
determined by inverting TL, 



b = —C 77,177-2 , 
C = — C 77-177,3 , 

P = -Cn 2 n 3 . 



,) = \-H- l \v), 



and turns out to be directed along the unit vector n 

1 - 

1 + eft"" 



I Poo/ = 



|n). 



(3.10) 



(3.11) 



Inserting these components in the expansion (3.6), one finds that the asymptotic density 
matrix is purely thermal, with a temperature given by the Unruh temperature: 



Poo = 



Trfe-ft^s] 



(3.12) 



Therefore, a two-level system which is uniformly accelerating in a vacuum scalar field is 
driven to a thermal state with temperature Tjj, irrespectively of its initial state: this 
thermalization phenomenon is the most obvious manifestation of the Unruh effect in the 
framework of open system dynamics. 

Nevertheless, further aspects of this phenomenon can be analyzed by studying the 
behaviour of the solution of (3.7) for finite times, that can be formally written as 



\p(t))=M(t)\p(0))+ 1-M(t) 



M(t) = e 



— v-^Ht 



(3.13) 



The matrix M.(t) is defined through the series expansion of the exponential function and 
therefore seems to involve arbitrary powers of TC. However, by definition, this 3x3 matrix 
obeys its cubic eigenvalue equation, so that powers of Tt higher then two can always be 
reduced to combinations of 7i 2 , H and 1, the unit 3x3 matrix. Then, a systematic use of 
this substitution allows to write: 



M(t) = 



e -^t Ai +2e -2(2 A+C)t 



A 2 cos Qt + A 3 



O 2 + 4C 2 

where the three constant 3x3 matrices Aj are explicitly given by: 

i-2(2A + c) n + n 2 , 
i 



sinOt 



(3.14) 



Ai = 
A 2 = 
A 3 = 2A 



(2A + C) 2 + ^ 



-2A(A + C) 1 + (2A + C) H - 
r 2 



n 2 



(3.15) 



-C(2A + C) 



1 + 



C(4A + C) - 



s 



n-cw 



n 



As expected, A4(t) contains exponentially decaying factors involving the real parts of the 
eigenvalues of 7i, modulated by oscillating terms in the effective frequency O. In other 
terms, an accelerating atom immersed in a vacuum scalar field is subjected to phenomena 
of decoherence and dissipation, all regulated by the Planckian factors appearing in the 
Kossakowski matrix (3.2), (3.3). 

These non-unitary effects can be studied by analyzing the time behaviour of suitable 
atom observables. Indeed, any physical property of the moving atom can be represented 
by an hermitian matrix O, that can be conveniently decomposed as 

3 

fj,=0 

The time behaviour of its corresponding mean value is then determined by that of the 
density matrix p(t): 

3 

(0(t))=Tr[Op(tj\ =O + Y t OiPi(t) . (3.17) 

i=i 

When the observable O represents itself an admissible atom state pi , the mean value (3.17) 
coincides with the probability Pi_»f(t) that the evolved atom density matrix p(t), initially 
in p(0) = pi, be found in such a state at time t. Using (3.14), this probability can be 
computed in general: 



1 - e~ PuU 
1 + e-P u " 



+ e-* At • ft) (pf • ft) 



+ e - 2 ( 2A + c )* ^ (p^ ■ pf ) - (p*i . fi) (p f • n) cos fit - n ■ (pi x p f ) sin ttt j | , 

(3.18) 

where, expanding the density matrices pi, pf as in (3.6), the notations (pi • pf) and (p; x pf) 
(and similarly with n) represent scalar and vector products of their corresponding coherence 
vectors. 

When pi = — n and p*f = n, the density matrices p;, pf represent the ground and 
excited states of the system hamiltonian Hs in (2.2) (see Appendix). In this case, the 
expression in (3.18) simplifies, 

giving the probability for a spontaneous transition of the atom from the ground state to its 
excited state. It is to this phenomenon of spontaneous excitation that one usually refers 
when discussing the Unruh effect; indeed, (3.19) vanishes as j3u — > oo, i.e. for an atom at 
rest. 

Although the behaviour of Vi^{t) in (3.18) and (3.19) is in principle experimentally 
observable through the use of suitable interferometric devices, in standard analysis of 
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the Unruh effect one limits the discussion to the spontaneous excitation rate i~i_>f, the 
probability per unit time of the transition i — > f , in the limit of an infinitely slow switching 
on and off of the atom-field interaction. In our formalism, its expression can be easily 
obtained by taking the time derivative of Vi->{(t) at t = 0; in the case of (3.19), one then 
finds 

/U f = - -=— , (3.20) 

which is the expected result for an interaction hamiltonian of the form (2.3). One should 
nevertheless remark that the possibility of a non-vanishing i~i_,f is just one of the many 
manifestations of the Unruh effect, that, as discussed above, actually involves phenomena 
of decoherence and dissipation; in this respect, the open system approach to the description 
of a moving atom dynamics appears to be much more physically comprehensive than the 
most traditional treatments. 



4. TWO ACCELERATING ATOMS AND ENTANGLEMENT 
ENHANCEMENT 

In the previous section we have seen that a uniformly accelerating two-level atom im- 
mersed in a scalar field in its vacuum state can be consistently described as an open system 
in weak interaction with a heat bath: in the comoving frame, the atom is seen evolving in 
time according to a master equation in Kossakowski-Lindblad form, that through decoher- 
ence effects drives its state towards a purely thermal equilibrium state, characterized by the 
Unruh temperature. When the system that is subjected to the uniform acceleration along 
the trajectory (2.1) is formed by two, non-interacting two-level atoms, one thus expects 
similar mixing-enhancing phenomena to occur, leading in particular to loss of the mutual 
quantum correlation (entanglement) that might have been present at the beginning. 

However, even though not directly coupled, the external vacuum field through which 
the two atoms move may provide an indirect interaction between them, and thus a means 
to entangle them. Indeed, entanglement generation through the action of an external heat 
bath has been shown to occur in certain circumstances; [18-23] it is therefore of physical 
interest to investigate the same issue in the case of accelerating atoms. 

We shall therefore start by considering a system composed by two, equal two-level 
atoms, that start moving along the trajectory (2.1) at proper time t = 0. Being inde- 
pendent, without direct mutual interaction, in the common comoving frame their internal 
dynamics can again be taken to be described by the generic hamiltonian (2.2). Then, the 
total two-system hamiltonian H$ is now the sum of the two terms: 

3 3 

H s = i4 1} + H { p , HP = ^J2 n ^® °o) , Hf ] = ^>(^ a t ) . (4.1) 

i=i i=i 

Similarly, being immersed in the same field $ M and within the weak-coupling hypothesis, 
the atom-field interaction hamiltonian can be most simply assumed to be the generalization 
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of that in (2.3): 

3 

H' = ® <ro) + {°o ® ^)] ■ ( 4 - 2 ) 

On the other hand, the field Hamiltonian H§, remains that of a collection of free, indepen- 
dent scalar fields. 

The derivation of the appropriate master equation describing the dynamics of the two 
atoms in the comoving frame proceeds as in the case of a single moving atom, discussed in 
Section 2. One starts from the Liouville-von Neumann equation (2.6) generating the time 
evolution of the state ptot(O) of the total system {atoms -f- external fields}, and then traces 
over the fields degrees of freedom, assuming a factorized initial state ptot(O) = p(0)<8> |0)(0|. 
In the weak-coupling limit, the two-atom system density matrix p(t) is seen evolving in 
time according to a quantum dynamical semigroup of completely positive maps, generated 
by an equation of Kossakowski-Lindblad form: 

^ = -i[H^p{t)]+C[p{t)]. (4.3) 

The unitary term depends on an effective hamiltonian H e ^ which is the sum of H$ in 
(4.1) and suitable Lamb contributions. In order to discuss them explicitly, we proceed as 
in the single atom case and adopt the simplifying condition (3.1), that results in diagonal 
field correlations (2.14). In this case, the effective hamiltonian consists of the sum of three 
pieces: H e s = H^g + + H^ 2 \ The first two represent single system contributions; 
they can be written exactly as in (4.1), with the frequency 10 replaced by the renormalized 
one O given in (3.5). The third term is a field-generated direct two-atom coupling term: 

3 

#£ 2) = 1 { I^M + /C (-^)] 6 v + M ) " ^M - rairij} a, <g> a 3 , (4.4) 

where /C(A) is the function introduced in (2.22). As explained at the end of Section 2, a 
suitable acceleration independent subtraction has implicitly been included in the definition 
(4.4) in order to make the contribution well defined. Further, recall that /C(A) can 

be split as in (2.23) into an acceleration dependent and an a = piece. Since, as observed 

(12) 

there, the acceleration dependent contribution to /C(A) is odd in A, one deduces that 
does not actually involves a: it is the same Lamb term that would have been generated in 
the case of a two-atom system at rest. Being interested in acceleration induced effects, and 
in particular in those related to entanglement creation, we shall not consider any further 
this inertially generated term, nor the single system contributions , and move to 
analyze the effects produced by the dissipative term C in (4.3). 

Since the interaction of the two moving atoms with the external scalar field is mediated 
by the same field operator $>(x), the Kossakowski matrix in C[p] involves the same field 
correlation functions (2.14) discussed in the case of a single atom system, and therefore can 
be expressed in terms of the single atom Kossakowski matrix in (3.2), (3.3). Explicitly, 
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one finds: [23] 



a 



(<Tj <g) CT ) /O ((7j <g) (7 ) - - 1 {<Ji<7 <g) (T ) , p| 

(cr ® aj) p(ao <S> - ^{(°"o <E> o-iCTj) , p| 

(<Tj <g> CT ) p(cr <S> <7») - (8) <Tj) , p| 

(cr <g> CTj)p(ai <g> cr ) - ^{(^i ® > p} 



(4.5) 



+ 
+ 
+ 



Describing the states of two two-level systems, the density matrix p{t) is now a 4 x 4 
matrix; in analogy with the decomposition (3.6) and recalling the normalization condition 
Tr[p(i)] = 1, we find convenient to decompose it as: 

]T 3 3 3 

p(t) = - a <g> cr + ^2 P^i 1 ) a o ® °i + E Pi0 ^ a i® a o+ E Pi i^ 0i ® °i ' ( 4-6 ) 

L i=l i=l i,J = l 

where the components poi(t), Pio(t), Pij(t) are all real. Substitution of this expansion in 
the master equation (4.3) allows deriving the corresponding evolution equations for the 
above components of p(t). As explained above, we shall ignore the hamiltonian piece in 
(4.3) since it can not give rise to acceleration induced entanglement. Further, we shall 
work in the regime of large acceleration, i. e. in the limit of flu small; while the conclusions 
concerning entanglement enhancement are actually independent from this choice, it will 
make more readable many explicit formulas. In fact, with this simplifying assumption, the 
Kossakowski matrix in (3.2) reduces to:" 1 " 

ciij = A Sij - iB e ijk n k , A= \ , B = ^- . (4.7) 

■iTTflu 47T 

A straightforward but lengthy calculation allows then derive from (4.3), and (4.5-7) 
the following result: 

3 

-AAp 0l (t) + B(l + 2T)n,-2B^2n kPlk (t) , (4.8a) 



dpoijt) 
dt 

dpiojt) 
dt 

9pij(t) 
dt 



k=i 

3 



- AAp l0 (t) + B(l + 2r) n,-2Bj2 n k Pki(t) 



(4- 



k=i 



= -4A[2p ij (t)+p ji (t)-r6 ij ] + 4B[m p 0j (t) + rij p i0 (t)] 



+ 2B[n i p j0 (t) + n j p 0i (t)] - 2B 5 2J ^ n k [p k0 (t) + p ok (t)] . (4.8c) 



k=i 



t Notice that the condition B < A, implicit in the original expressions for A and B in 
(3.3), needs now to be formally imposed in order to maintain the positivity of a^-; it is 
physically justified by the assumption of a small flu- 
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In these formulas, the quantity r = X^=i Pa represents the trace of p^j ; it is a constant of 
motion, as easily seen by taking the trace of both sides of (4.8c). Further, the value of r can 
not be chosen arbitrarily, since it has to comply with the requirement of positivity of the 
initial density matrix p(0); indeed, using the decomposition (4.6), one finds: — 3 < r < 1. 

The system of first order differential equations in (4.8) naturally splits into two inde- 
pendent sets, involving the symmetric, p( i) = Poi+Pio, P(ij) = Pij+Pji, and antisymmetric, 
P[0i] = Poi — PiQi P[ij] = Pij — Pji-, variables. Although both sets of equation can be exactly 
integrated, the form of the explicit solutions looks cumbersome and is not very inspiring. 
Nevertheless, by just examining the structure of the two sets of differential equations, one 
can conclude that the antisymmetric variables involve exponentially decaying factors, so 
that they vanish for large times, while the symmetric variables approach in the same limit 
a non-vanishing asymptotic value. As a consequence, the evolution equations (4.8) admit 
an equilibrium state p, whose explicit form can be obtained by inverting the appropriate 
coefficient tensors multiplying the variables poi(t), Pio(t), Pij(t) in the r.h.s. of (4.8). 

The equilibrium density matrix p can be expanded as in (4.6): its components are 
given by: 

POi = PiO = 3 + R2 ( T + 3 ) Ui ' 

(4.9) 

1 



PlJ ~ 3 + tf 2 



{t-R 2 ) Sij +R 2 (t + 3) 



where R = B/A is the ratio of the two constants appearing in the Kossakowski matrix in 
(4.7), whose positivity implies: < R < 1. As expected, the antisymmetric components 
P[ j] and /3[jj] are zero, while the only dependence on the initial state is through the constant 
t = J2iPu(0)^ It is remarkable that in general this equilibrium state turns out to be 
entangled. 

In the case of two, two-level systems, a measure of the entanglement content of any 
state p is provided by the concurrence C[p\. [41-43] Indeed, it has been shown that C is a 
monotonically increasing function of the entanglement of formation; [44] its value ranges 
from zero, for separable states, to one, for fully entangled states, like the Bell states. In 
order to compute the concurrence of any 4x4 density matrix p representing the state of two 
atoms, one starts from the auxiliary matrix p = (02 <S> cr 2 ) p T (02 ® 02), where T indicates 
transposition. Although not necessarily hermitian, the matrix pp has real, non-negative 
eigenvalues, whose square roots A M , p, = 1, 2, 3, 4, can be ordered in decreasing order: Ai > 
A2 > A3 > A4. The concurrence of pis then defined to be: C[p] = max{Ai — A2 — A3 — A4, 0}. 

In the case of the asymptotic state p in (4.9), one finds that the concurrence is indeed 
non-vanishing provided 

5R 2 -3 , A . 

r < -j^rw • (410) 



t One can check that the form (4.9) for the equilibrium state remains unchanged even 
considering the more general evolution generated by the Kossakowski matrix (3.2) instead 
of the simplified version in (4.7). The considerations about entanglement production below 
are therefore valid in general and not only in the case of large accelerations. 
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and that, in this case, its explicit expression is given by: 



C[ P ] 



(3 - g) 
2(3 + R 2 ) 



~5R 2 -3 
3-R 2 



(4.11) 



The concurrence is therefore a linearly decreasing function of r, starting from its maximum 
C[p] = 1 for r = — 3 and reaching zero at r = (5R 2 — 3)/(3 — R 2 ); notice that this ratio 
is an admissible value for r, since it is always within the interval [—1,1] for the allowed 
values of R. 

This result is remarkable, since it implies that the dynamics in (4.8) can generate 
entanglement: one prepares the two atoms in a separable state at t = 0; then, provided 
the condition (4.10) is satisfied, their long time equilibrium state will turn out to be 
entangled. The simplest example of a separable state is provided by the direct product of 
pure states: 



P(0) 



1 , 

- [ 1 + m • a 

At 



(4.12) 



where n and m are two unit vectors. In this case, one easily finds that r = n • m, so that, 
recalling (4.11), the asymptotic entanglement is maximized when n and rh are collinear 
and pointing in opposite directions. Explicitly, one finds: 



C[p} = 



2R 2 



3 + R 2 ' 



(4.13) 



which reaches its maximum value, C[p] = 1/2, for R = 1. 

The phenomenon of entanglement production occurs also in cases when the initial 
state p(0) already has a non- vanishing concurrence. As observed before, C[p] reaches its 
maximum when r = —3; in this case, the equilibrium state p coincides with the totally 
entangled singlet state: 

1 [ 3 1 

P- = ^ oo <S> cr - ^o-j <S> cfi , (4-14) 

L i=i 

a fixed point of the dynamical equations (4.8), as easily seen by direct inspection. Then, 
let us consider the following initial state, 



p(0) = (1 - e)p- + -<j <g> cr , 



(4.15) 



which interpolates between p_ and the totally mixed (separable) state; for e < 2/3, it is 
entangled, with C[p(0)] = 1 — 3e/2. In passing from this initial state to its corresponding 
equilibrium state p as t — > oo, the corresponding increase in concurrence, and thus of 
entanglement, can be easily computed: 



C[p] - C[p(0)} = 



3R 2 e 
3 + R 2 



(4.16) 
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which is indeed non vanishing. Further, notice that for the state (4.15), r = —3(1 — e), so 
that it can be taken to be very close to its lower limit —3; unfortunately, the entanglement 
production (4.16) becomes vanishingly small as p(0) approaches p_. In other terms, the 
maximally entangled state p_ can never be asymptotically reached, and thus the maximum 
of concurrence obtained, unless one already starts with it at t = 0. 



5. DISCUSSION 

Moving detectors, modelled as simple two-level atoms, immersed in external vacuum 
fields and following an uniformly accelerating trajectory are seen to posses a non-vanishing 
probability of spontaneous excitation, reproducing a thermal spectrum: this phenomenon 
is usually referred to as the Unruh effect. 

It turns out that the dynamics of these accelerated atoms can be consistently assimi- 
lated to that of subsystems in interaction with an external environment, i.e. the so-called 
open quantum systems. General techniques and results obtained in the analysis of the latter 
can then be fruitfully applied to the study of the former, enabling the discussion of physical 
aspects of the Unruh effect that might not be easily identified in the standard treatments. 
As shown in the previous sections, the open system paradigm is particularly suitable for 
analyzing on one hand issues connected to the appearance of decoherence effects, and on 
the other hand questions related to the phenomenon of entanglement enhancement. 

The starting point of our analysis has been the derivation of the appropriate master 
equation that generates the time evolution of the states of an accelerating system, in its 
comoving frame. In the framework of a physically justified weak-coupling hypothesis, the 
system subdynamics takes the form of a semigroup of completely positive maps, trans- 
forming density matrices into density matrices, while preserving their normalization and 
positivity. 

In the case of a single two-level atom, the master equation has been explicitly inte- 
grated; this has allowed discussing in detail the mixing-enhancing properties of the asso- 
ciated finite time evolution, through the analysis of the behaviour of appropriate atom 
observables. In particular, one finds that the moving atom is subjected to dissipative ef- 
fects, that asymptotically drive its density matrix to an equilibrium state, with a purely 
thermal spectrum. t 

When the accelerating system is composed by two, independent atoms, its physical 
characterization naturally involves the analysis of their mutual quantum correlations. Be- 
cause of the just mentioned decohering and mixing-enhancing phenomena, one would be 
led to regard the Unruh effect as counteracting entanglement enhancement. 

Nevertheless, quite in general the presence of an external environment (in the specific 
case of the external fields) can provide indirect interaction between the two otherwise 

t Note that this result is different from (but complementary to) the "thermalization 
theorem" discussed in the usual treatments of the Unruh effect: there, it is the field 
"vacuum" state that appears to be purely thermal to an accelerating observer. 
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totally decoupled two-level systems, thus a means to correlate them. This picture has 
indeed been confirmed by the analysis of specific models; in particular, [23] entanglement 
can be created by the action of an external bath through a purely noisy mechanism during 
the memoryless, Markovian regime, when the corresponding subdynamics is generated by 
equation of the form (2.11), (2.12). 

The master equation describing the dynamics of two accelerating atoms discussed in 
Section 4 is precisely of the form identified in Ref.[23] as generating initial entanglement. It 
is remarkable that the unavoidable decoherence that subsequently builds up is not sufficient 
to counteract this quantum correlations enhancing effect. The fate of the entanglement 
production by acceleration can in fact be discussed by analyzing the entanglement content 
of the final two-atom equilibrium state, through the evaluation of its concurrence: quite 
in general, one finds an asymptotic non-vanishing concurrence even for totally separable 
initial states. 

Although obtained in an idealized situation, this result offers new possibilities for 
an actual direct verification of the Unruh effect. So far all efforts have been devoted to 
the analysis of possible experimental settings that could allow the measure of the tiny 
spontaneous excitation thermal rate induced in single accelerating systems. Alternatively, 
using suitable devices, one can instead try to detect the quantum correlation enhancement 
that is generated when the accelerating system is formed by two, initially un-entangled 
atoms; in view of the high accuracy and sophistication reached by present quantum optics 
experiments, this possibility might actually be realized in the future. 

As a concluding comment, let us remark that the presented open system treatment of 
the Unruh effect is not limited to the analysis of the simple setting of atoms interacting 
with free, scalar fields. Extension to the case of higher spin fields is straightforward, while 
the analysis of situations involving self-interacting fields would require the use algebraic 
quantum field theory; this might not be such a formidable task as it looks, since the 
quantum theory of open systems is amenable to a rigorous, algebraic formulation. Finally, 
it is known that the Unruh effect has many similarities and analogies with the phenomenon 
of particle creation in curved spacetimes. In this regard, we expect the open system 
paradigm to be applicable also to those cases, possibly providing new insight in the physical 
interpretation of the corresponding effects. 
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APPENDIX 

Using the rigorous results of Refs.[16, 17], we shall here present the derivation of 
the master equation generating the reduced dynamics of a two-level system in interaction 
with a generic external environment, in the weak coupling limit. We shall generalize the 
treatment of Ref.[45] and take the generic hamiltonian (2.2), 

H s = ^n-a, |n| = l, (Al) 

as system hamiltonian; instead, the most general term that is linear in both system and 
environment variables, 

3 

ti=0 

will be taken to represent the interaction hamiltonian. The environmental dynamical vari- 
ables are left unspecified (they coincide with the field in the case of an accelerating 
system), as the corresponding hamiltonian ifg, that nevertheless is assumed to generate 
time translations: 

e lHBt B^e~ lHBt = B„(t) . (A3) 

The time evolution of the density matrix p tot representing the state of the complete 
system is then generated by the total hamiltonian 

H = H + H' , H = H s <g> 1 + 1 <g> H B , (A4) 

through the standard unitary evolution, starting with an initial state taken in factorized 
form: ptot(O) = p(0) <8> ps', p(0) is the 2x2 density matrix describing the subsystem state, 
while pb is the analogous one for the environment, assumed stationary: [Hb,Pb] = (in 
the case of moving atoms, ps = |0)(0|, with |0) the Minkowski vacuum). 

Correspondingly, the subdynamics describing the evolution of the subsystem alone 
is obtained with a trace operation over the environment degrees of freedom. As briefly 
described in Section 2, and rigorously proven in Refs.[16, 17], in the limit of weak coupling 
between subsystem and environment, the reduced density matrix p(t) = Tr[p to t(t)] is found 
obeying the following evolution equation: 

^ = -iL Hs [p{t)\ + K) [p(t)\ , L Hs [p\ = [H S , p] , ( A5) 

where 

£"[•] = - lim 1 / ds U{-s) JCU{s)[-} , U{s) = e ~ isLH s , (A6) 

T^oo T Jo 

JC[p] = J™ dt Tr([e* Hot H> e-* Hot , [H' , p®p B ]]) ■ (A7) 

The ergodic mean over the system dynamics in (A6) is necessary in order to obtain a 
completely positive subdynamics, and it is physically justified as an averaging over the 



and 
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fast microscopic system oscillations; indeed, the weak coupling limit procedure involves 
a suitable rescaling of the time variable, so that the equation (A 5) actually generates a 
"coarse grained" subdynamics with respect to the free system motion. With the definitions 
(A.l) and (A. 2), one can more explicitly write: 

/C«[p]= lim - ds dt <T li (t + s)p<T v (s)(B v B li (t)) + <T v (s)p<T IM (t + s) (B^Bu) 
t^oo 1 Jq Jo l 

+ s) a v {s) p (B^t) B v ) - p a u {s) a^t + s) {B v , 

(A8) 



where 



(B^t)B v ) = Tr[B IM {t)B vPB ] , 



(A9) 



are the environment correlations. 

In order to proceed further, it is convenient to introduce the two projectors operators 



P ± = i(i±*.*) 



(A 10) 



they represent the density matrices of the two eigenstates of the system hamiltonian (A.l), 
with eigenvalues ±u/2. One can then use the auxiliary matrices a^P, £ = 0, +, — , explicitly 
defined by 

4°) = P+ a M P+ + P. a M P. , <t& = P± a M P T , (All) 

and represent the system free time evolution in terms of the following spectral decompo- 
sition: 

a^t) = e lHst a ll e- lHst = ^ e* ut a® . (A12) 

€=o,± 



This allows performing explicitly the limit in (A8) and expressing the result in terms of 
the following Fourier and Hilbert transform of the environment correlations: 



/oo 
dte^ (0|fl M (f)5„|0) , 
-oo 



(A13) 



and 



/•OO />oo 

= / dt e* ut (0|S M (t)^|0) - / dt e~* ut (0\B^B„(t)\0) ; (A14) 
Jo Jo 

notice that the first 4x4 matrix is hermitian, [ajLD (w)]^ = while the second is 

antihermitian, [(3$ (a;)] ^ = —f3$(u). Explicitly, one finds: 



£=0,± M,^=0 



P , 



(A15) 
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This expression can be further simplified by expanding the auxiliary matrices in terms 
of Pauli matrices: 



with 



/(o) 



(A16) 



(A.17) 



The entire master equation in (A5) can then be rewritten in standard Kossakowski- 
Lindblad form: 

Mf) =-i[H eBi p(t)] + C\p(t)] , (A18) 



dt 



where 



H e ft 



1 3 

*,J = 1 
1 3 

- ^ [w n t + hi] , 



(Al9a) 
(A 196) 



i=l 



while the Kossakowski matrix and hamiltonian vector 6j can be expressed as: 



ciij 



E E 



(A20a) 



£=0,± fc,«=l 
3 



(0) 
0j 



(0) 



«(0) _ o(0) 



n 7 - rii 



+ ^2 

j,k,l,m=l 



E 

€=o,± 



(A206) 



Since [V^]* = ' the matrix a^- in (A. 20a) results manifestly hermitian and posi- 

tive, being the combination of Fourier transform of correlation functions. [14] Further, by 
choosing for the environment variables the fields $ M , the master equation (A 18 — 20) 
reduces to that discussed in the text. 

Notice that in general the dissipative piece in (A. 19a) is a function of nine, real 
parameters, the independent entries of the matrix ciijJ However, because of the structure 
given in (A. 20a), the Kossakowski matrix obtained through a weak coupling procedure 
appears to depend on a lesser number of free parameters. This fact has been already 
noticed in Ref.[45], where it is further observed that the most general master equation 



t These are not actually all independent: since the matrix ciij is positive, they need to 
satisfy certain inequalities; see Ref.[28] for details. 
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obtainable in the weak coupling limit appears to coincide with the old Bloch equation, 
describing the dissipative motion of a spin in a constant magnetic field. 

In the light of the previous derivation, this conclusion looks however too restric- 
tive: by starting with the most general system hamiltonian (Al), instead of the one 
with n = (0, 0, 1) as adopted in Ref.[45], one is able to obtain a master equation in the 
form (A 18 — 20), certainly more general than the Bloch equation. This observation might 
have interesting applications in the study of specific open system models. 
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